LAW OF THE ABSORPTION TIME OF POSITIVE SELF-SIMILAR 

MARKOV PROCESSES 



P. PATIE 



Abstract. Let X be a positive self-similar Markov process with as an absorbing state. 
The purpose of this paper is to describe the law of the absorption time, say To, which might 
occurs continuously or by a jump. We start by showing that the distribution function of To 
can be expressed in terms of an increasing invariant function for a specific transient Ornstein- 
Uhlenbeck process associated to X. Furthermore, specializing on the spectrally negative case, 
we suggest an original methodology to get a power series or an integral representation of this 
invariant function. Then, by means of probabilistic arguments, we deduce some interesting 
analytical properties satisfied by these functions, which include, for instance, several types of 
hypergeometric functions. We end the paper by detailing some known and new examples. 
In particular, we offer an alternative proof of the recent result obtained by Bernyk et al. 2 
regarding the law of the maximum of regular spectrally positive stable processes. 



1. Introduction 

Let X = ({Xt)t>o, (Q x )x>o) be a self-similar Hunt process with values in [0, oo). It means 
that X is a right-continuous strong Markov process with quasi-left continuous trajectories and 
there exists a > such that X enjoys the following self-similarity property: for each c > and 
x > 0, 



the law of the process (c 1 X c a t ) t > , under Q x , is 



tx te- 



at is called the index of self-similarity. This class of processes has been introduced and studied 
by Lamperti [23]. During the last two decades, they have been intensively studied, to name just 
a few, we mention the works [15], [15], [7J, [11], [3D], [TTJ and [37]. Furthermore, we recall that 
Lamperti's main motivation was the fact that self-similar Markov processes arise as the limit 
of properly normalized stochastic processes. It is then not very surprising that these processes 
appear in different settings, such as, for instance, coagulation- fragmentation processes [5], and 
continuous state branching processes with immigration [35]. We also mention that, recently, in 
the SLE context, Alberts and Sheffield pQ describe a measure-valued function supported on the 
intersection of a chordal SLE(k) curve with R, 4 < k < 8, in terms of the law of the absorption 
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time Tq of some Bessel processes. Next, Lamperti [21] states the following zero-one laws: if 

T = inf{s > 0;X S = 0} 

then, for all x > 0, either Q x (To < oo) = 1 or Q x (To < oo) = 0. In this paper, our purpose 
is to describe the distribution function of To in the former case. Unfortunately, beside some 
isolated cases the distribution of To is not attainable. We mention the papers [2], [19] and [35] 
where such examples can be found and refer to the survey paper [8] for a description of these 
cases. Besides, two notable exceptions might be worth mentioning: when X is a Bessel process 
of negative index and when X is a regular spectrally negative Levy process killed upon entering 
the negative half-line. In the former case, several proofs can be found in the literature, see for 
instance the excellent monograph of Yor [45J and the more recent survey papers of Matsumoto 
and Yor [26] and [2Z]. However, most of the proofs rely on the knowledge of the semigroup of 
Bessel processes. For the second case, Bernyk et al. [2j derive a representation of the distribution 
of To by inverting, in a non trivial way, the known expression of the Wiener- Hopf factorization of 
regular stable one-sided Levy processes. Our approach will differ from these two cases. Indeed, 
in general, we do not have access neither to the semigroup of X nor to the Laplace transform 
of T). We shall first show that the distribution function of To can be described as an invariant 
function for a specific Ornstein-Uhlenbeck process associated to X. Then, specializing on the 
case when X has only negative jumps, we suggest a methodology, relying merely on probabilistic 
arguments, to represent this invariant function as a power series. It turns out that finding the 
distribution function of To boils down to compute merely an exponential moment for the first 
passage time above for X. 

We now take X to be realized as the coordinate process X t {oj) = u>(t) on the Skohorod space 
of cadlag functions from [0,oo[ to itself. We assume that is an absorbing state for X, so 
that each of the laws Q x governing X is carried by {u> E £1 : oj(t) = 0, Vi > To (a;)}. The natural 
filtration on Q is (Ft)^, and J 7 ^ = Ut>o^ r *- Lamperti proved that for each fixed a > 0, there 
is a bijective correspondence between [0, oo)-valued self-similar Markov processes with index a 
and real- valued Levy processes, i.e. processes with stationary and independent increments. More 
specifically, by introducing the additive functional 

E t = inf{s > 0; A s = [ X~ a dr >t}, t>0, 
Jo 

Lamperti [21] showed that the process £ = (&)t>o, defined by 
(1-1) it = log(X Et ), t>0, 

is a Levy process. Observe that the additive functional (A t )t>o is then determined by 

A t = inf{s > 0; S s = [ e a ^ r dr >t}, t> 0. 
Jo 

We write P for the law of the Levy process with £o = 0- Moreover, let be the characteristic 
exponent of £. It is well known, see [3] and [22] for background on Levy processes, that it admits 
the following Levy-Khintchine representation, for any u £ zK, the imaginary line, 

(1.2) = bu + -u 2 + / (e ur - 1 - url { \ r \ <l} )u{dr) - q 

2 Jr 

where the parameters q > 0,b £ M.,a > and the measure v which satisfies the integrability 
condition L(l Ar 2 ) v{dr) < +oo uniquely determined the law of £i. We shall refer to £ (resp. \&) 
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as the underlying Levy process (resp. characteristic exponent) of X. Furthermore, in the case 
Q x (Tq < oo) = 1, Lamperti [24] explained that, either q > and the absorption time occurs by 
a jump, i.e. 

Q x (X To _ > 0, T < oo) = 1, Vx > 0, 

or £ drifts to — oo and Tq occurs continuously, i.e. 

Qx(X To _ = 0, T < oo) = 1, Vx > 0. 

By means of the strong law of large numbers for Levy processes we gather these two possibilities 
in the following. 

HO. Either q > or lim^oo ^ = — oo a.s. and q = 0. 

We emphasize that throughout the paper this condition does not hold necessarily, that is we 
may also consider the case when Q x (To < oo) = 0,x > 0. 

The remaining part of the paper is organized as follows. In the next Section, we show that 
the distribution function of To can be expressed in terms of an increasing invariant function 
for a transient Ornstein-Uhlenbeck process associated to X. We also recall several connections 
between the law of To and some other important positive random variables. In Section [3l 
specializing on the case when X has only negative jumps, we derive a power series and an 
integral representation of the distribution of To. The proof and several consequences of this 
result are also presented in this part. Finally, in the last Section, we detail some known and 
new examples. We also mention that some of the results stated in Theorem 13.21 below were 
announced without proofs in the note [38] . 

2. A general result and some interesting connections 

Throughout this Section, we assume that HO holds. We define the Ornstein-Uhlenbeck process 
U = (Ut)t>o associated to X by 

(2.1) Ut = e &t X {1 _ e - t] , t >0, 
where a = a . Next, we introduce the stopping time 

(2.2) Ho = inf{s > 0; U s = 0} 
and the function 

1 _ p[x) = Q X (H < oo), x > 0. 
We are now ready to state the following. 
Theorem 2.1. For any x > and t > 0, we have 

(2.3) P(xt~ & ) = Q x (T >t) 
and P is increasing on M + with linx^oo P(x) = 1 and P(0) = 0. 

Proof. First, a simple time change yields the following identity in distribution 

Fo = -log(l-T Al). 
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Thus, we deduce that 

l-P(x) = Q z (#o<oo) 
= Q X (T <1). 

Then, invoking the self-similarity property of X we obtain the identity 

Qx(T >t) =® xt - & {T > 1) 
from which we deduce the identity (|2.3p and the properties stated on P. □ 



We proceed by providing an heuristic analytical argument explaining the identity (|2.3p . We 
need to introduce the notation 

Q(t,x) = Q x (T >t), x,t>0, 
and to state the following result which is found in [24|, Theorem 6.1]. 

Proposition 2.2. Let f : R + — ► M be such that f(x),xf'(x) and x 2 f"(x) are continuous func- 
tions on R + , then f belongs to the domain, O(L), of the characteristic operator L of (X, Q) 
which is given, for any x > 0, by 

(2.4) L/(s) = x" Q (|x 2 /" (x) + 

+ f{rx) - f{x) - xf'(x) log(r)I { | log(r) |< 1} i/(dlog(r)) - qfix) \ . 

Moreover, if HO holds, then f G B(L) if and only if f(0) = 0. 

Note that the process U defined in (|2.ip is a transient Hunt process on [0, oo) and its infini- 
tesimal generator, denoted by A, has the following form, see e.g. Carmona et al. [14J, 

Af(x)=-Lf(x)+axf'(x), x>0, 

for a function / G D(L). Next, let us assume that the mapping it, x) i— > Q(t,x) defined on 
M + x M + is smooth enough and thus can be expressed as the solution to the following initial 
boundary value problem 

{^Qit,x) =~LQ(t,x) on (t,x) G (0,oo) x (0, oo), 
Q(t,0) =0,t>0, 
Q(0,x) =1, x>0. 

As above, from the self-similarity property of X, one gets that 

Q(t, x) = Q(l, xt~ & ) = Qix~ a t, 1), t, x > 0. 

Thus, one may set 

P{xt~ & ) = Q(l,xt- & ), x,t>0. 

Next, observing that 

%-Q(t,x) = -axt- & - 1 P'(xt- & ) and ^-QCt,x) = t- & P'ixt~ & ), 
at ox 



one gets 



d d 
—Qit,x) = -axt~ l -^Qit,x) 



and the problem (|2.5p is equivalent to 

d_ 

i:)x 



LQ(t, x) + obi -1 x) = on (i, x) G (0, oo) x (0, do), 



Q(t,0) = 0, t > 0, 
Q(0,a;) = 1, x > 0. 

Note that the linear operator L acts on the space variable only. Hence, one may fix t = 1 and 
solve 

'AP(i) = onxe(0,oo), 
P(0) = o, 
[P(oo) = l. 

Therefore, we have transformed a parabolic integro-differential initial-boundary value problem 
into an integro-differential equation with boundary conditions. Actually, such a transformation is 
well-known for some specific partial differential equations. Indeed, by using a linear combination 
of the one-parameter groups of scaling symmetries of the heat equation, one obtains as reduced 
equation the Weber's differential equation, see [30] for more details. The solutions in such a case 
are called the scale-invariant solutions. 

2.1. Connection between To and some random variables. In this part, we establish or 
recall some connections between the law of the positive random variable To and the law of some 
interesting random variables. To this end, we introduce the following notation, which might be 
justified again by means of the self-similarity property of X, 

(2.6) S(tx~ a ) = P(xr & ) = Q X .(T > t), t,x>0, 

where we recall that a = a -1 . 

2.1.1. Law of the maximum of stable Levy processes. Let Z be an a-stable Levy process, with 
< a < 2 and we exclude the case when Z is a subordinator. Then, let us denote by X 
the process Z killed upon entering the negative half-line. X is then a positive self-similar 
Feller process and the characteristic triplet of the underlying Levy process has been computed 
by Caballero and Chaumont |12j . see also [13] and [23], for an expression of its characteristic 
exponent in terms of the Pochhammer symbol. Next, we denote by Z the dual of Z, i.e. Z = —Z 
which is also an a-stable Levy process. Then, by means of the translation invariance of Levy 
processes, we deduce readily the following identities 

Qx (T <t) = ¥ x ( inf Z s < 



max Z s > x 

0<s<t 



which can be written as follows 



(2.7) P max Z s > x ) = P(xt~ a ), x,t>0. 

\0<s<t J 

As a Levy process, it is well known that the law of the maximum process of stable processes can 
be described in terms of the Wiener-Hopf factors which corresponds to its Laplace transform 
in time and space. The methodology we will develop requires merely to inverse the Laplace 
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transform in space to obtain the law of the maximum process. We refer to Section [4.21 for more 
details on this case. 



2.1.2. Exponential functional of Levy processes. Let us recall the notation 

E t = f e a ^ 3 ds, t > 0. 
J o 

It is plain, from the Lamperti transformation (jl.ip . that, in the case q = and E[£i] < 0, one 
gets the identity in distribution 

Moreover, when q > 0, one has 

(T ,Q,)=V£ eg ,P) 

where e q is an independent exponential random variable of parameter q. Hence, we have for any 
t > and q > 0, 

(2.8) P(S e? >t) = S(t) 

where we understand that eo = oo. As pointed out by several authors, see Carmona et al. |14j . 
Rivero [3D], and, Maulik and Zwart [28], the study of the exponential functional is also motivated 
by its connection to some interesting random affine equations. Indeed, from the strong Markov 
property for Levy processes, which entails that for any finite stopping time T in the filtration 
C^Et) t — 0)) the process (£t+T — £t> t > 0) is independent of Ty, t and has the same distribution 
than £, we deduce readily that the random variable E e<3 is solution to the random affine equation 

(2.9) S eg ( £ } fe-t-ds + e-^ 

Jo 

where, on the right-hand side, is taken independent of (f^ e~^ s ds, e^ 1 ). This type of random 
equations has been deeply studied by Kesten [21] and Goldie [20]. Relying on Kesten's results, 
Rivero [40] shows that there exists a constant C > such that one has the following asymptotic 
behavior 

(2.10) S(t) ~ Cr M as t — > oo 

(f(x) ~ g(x) as x — > a means that lim x ^ a = 1 for any a £ [0, oo]) whenever the Levy 
process satisfies this set of conditions 

(1) £ is not arithmetic (i.e. its state space is not a subgroup of rZ for some r > 0), 

(2) there exists 6 > such that E[e e ^I { i <eq} ] = 1, 

(3) J x>1 xe 6x v{dx) < oo. 

We refer to Theorem l3.2l and Proposition l3.5l for several characterizations of the Kesten's constant 
C in the case £ is a spectrally negative Levy process. We end up this part by recalling that 
the exponential functional of Levy processes appears to be a key object in several fields. In 
particular, in mathematical finance, the law of S eq allows us to derive the Laplace transform in 
time of the price of the so-called Asian options, see [18] for the Brownian motion case and [33] 
for more general Levy processes. 
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2.1.3. Entrance law and stationary measures. Let us denote by X the dual of X with respect 
to the reference measure x a ~ 1 dx. Bertoin and Yor [6] showed that J is a positive self-similar 
Markov process with underlying Levy process £ = — £, the dual of ^ with respect to the Lebesgue 
measure. Moreover, assuming that £ is not arithmetic and oo < E[£i] < 0, they show that X 
admits an entrance law at 0, i.e. the family of probability measure (Q x )x>o converges, as x — > 0, 
in the sense of finite-dimensional distribution, to a probability measure, denoted by Qo+- In 
particular, they obtain the following identity, for any bounded borelian function /, 

Eo+ [/(*?)] = -^^^[/(S- 1 )^- 1 ]. 
a\E[^\\ 

Hence, we obtain for any x > 

Q 0+ ( Xl edx) = — ±— S(dx~ a ). 
\E[^i]\x 

We also mention that in [34J, it is shown that the measure Q +(^i € dx) is the stationary 
measure of the Ornstein-Uhlenbeck process defined by 

= e"'X et , t > 0. 

Finally, we define, under P, the process Y = (Yt)t>o, by 

Y t = e aft (y + J e'^ds) ,t>0. 

This generalization of the Ornstein-Uhlenbeck process is a specific instance of the continuous 
analogue of random recurrence equations, as shown by de Haan and Karandikar [16J . Carmona et 
al. [13] showed that Y is a homogeneous Markov process with respect to the filtration generated 
by £. Moreover, they showed, from the stationarity and the independency of the increments of 
£, that, for any fixed t > 0, 

Y t { ^Y e a ^+ [ e a ^ds. 
Jo 

Then, if E[fr] < 0, they deduced that Y t y -4Y, oo as t — > oo. Hence if Vo — Sqo, the process Y 
is stationary with stationary measure the law of X^. 



3. The spectrally negative case 

In this Section we consider the case when £ is a spectrally negative Levy process, i.e. v(0, oo) = 
in (jl.2p and we exclude the cases when £ is degenerate, that is when £ is the negative of a 
subordinator or a pure drift process. In this setting, it is well-known that the characteristic 
exponent admits an analytical continuation to the lower half-plane. Thus, we write simply 
ip(u) = fy(—iu),u > 0, when q = and ip{u) = ip{u) — q otherwise. It means that, for any u > 0, 
we have 

(3.1) i>(u) = bu + -u 2 + (e ur - 1 - url { \ r \ <l} )v{dr). 

^ J —oo 

Let us now recall some basic properties of the Laplace exponent ip, which can be found for 
instance in Bingham [10J. First, it is plain that lim M ^ 00 ^(ti) = +oo and by monotone con- 
vergence, one gets E[£\] = b + f_ ru{dr) G [—00,00). Thus, by means of the strong law of 
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large numbers for Levy processes, the condition lim^oo = — oo a.s. in HO is equivalent to the 
requirement E[£i] < 0. 

Moreover, the asymptotic behavior of if) depends on its coefficients. Indeed, if a = and 
J ^ I A |r|i/(dr) < oo, that is the Levy process £ has paths of bounded variation, then, as we 

have excluded the degenerate cases, we have, see [HI Corollary VII. 5], writing b = b — rv(dr), 
(3.2) lim = b > 0. 



+ 00 



For the other cases, we have 



(3.3) lim = +00 . 

u— >oo ii 

Differentiating ip again, one observes that if) is convex. Note that is always a root of the 
equation if>{u) = 0. However, in the case E[£i] < 0, this equation admits another positive root, 
which we denote by 6. This yields the so-called Cramer condition 

E[e^] = 1. 

Moreover, for any E[£i] G [— oo,oo), the function u i— ► ip{u) i s continuous and increasing on 
[max(#, 0), oo). Thus, it has a well-defined inverse function eft : [0, oo) — > [max(#, 0), oo) which is 
also continuous and increasing. Finally, we set the following notation, for any u > 0, 

ip$(u) = ip(u + 9), 

^<t>( q )(u) = ip{u + (f>(q)). 

Recalling that tp(9) = and observing that Tp((f)(q)) = 0, the mappings ipe^Mq) are plainly 
Laplace exponents of conservative Levy processes. We also point out that ip' e (0 + ) = ip'(6 + ) > 
and ?/>^ 9 )(0 + ) = ip'{4>{q)) = ^rj^y > 0. We denote the totality of all functions if) of the form (|3,ip 

with the additional condition that ip'(0 + ) > by CKL+. Note, from the stability of infinitely 
divisible distributions under convolution and convolution powers to positive real numbers, that 
£/C+ forms a convex cone in the space of positive valued functions defined on [0, oo). 
We proceed by introducing a few further notation. First, let if) E CK, + and set oq = 1 and for 
any n = 1, 2, . . ., 



a n (ip;a) = f if)(ak) j 
In [37], the author introduces the following power series 

oo 

(3.4) I^z) = Y J ^a)z n 

n=0 

and shows by means of classical criteria that the mapping z i— > T^(z) is an entire function. 
We refer to [37] for interesting analytical properties enjoyed by these power series and also for 
connections with well-know special functions, such as, for instance, the modified Bessel functions 
and several generalizations of the Mittag-Leffler function. The author also proves that for any 
r > 0, the mapping x i— * 2^(x a ) is an 1-eigenfunction for X and 

E r e -^l = ^fl, 0<x<a, 



where T a = inf{s > 0; X s = a}. Next, let G K be a Gamma random variable independent of 

—y k>— 1 

X, with parameter n > 0. Its density is given by g{dy) = 6 r ^ — c?y,y > 0, with T the Euler 
gamma function. Then, in [341 . the author suggested the following generalization 



(3.5) T^k-z) = E[lip(G K z)} 

1 r°° 

By means of the integral representation of the Gamma function T(p) = e~ s s p ~ 1 ds, 0\e(p) > 0, 
and an argument of dominated convergence, one obtains the following power series representation 

1 x 

(3.6) 1^{p;z) = — -^2a n {^;a)T(p + n)z n 

W> n =0 

which is easily seen to be valid for any \z\ < R = lim u _ +0O t^tl > q. Moreover, for any \z\ < R, 
the mapping p Z^(p; z) is a meromorphic function defined for all complex numbers p except 
at the poles of the Gamma function, that is at the points p = 0,-1,... However, they are 
removable singularities. Indeed, for any \z\ < R and any integer N E N, one has, by means of 
the recurrence relation T(z + 1) = zT(z), 2^,(0; z) = 1 and 

n=0 ' 

Thus, by uniqueness of the analytical extension for any \z\ < R, T^(p;z) is an entire function 
in p. Note also that for p = 0, —1, . . ., as a polynomial, I^(—p; z) is an entire function in z. We 
point out that when ip is the Laplace exponent of a standard Brownian motion with a positive 
drift, (I^p(—N; e l7r .z))jVeN is expressed in terms of the sequence of Laguerre polynomials. In the 
following, we summarize the above results and provide an analytical continuation of 1^{p; z) in 
the case R = ab, that is when £ is with paths of bounded variation. 

Proposition 3.1. (1) IflvaXu-^^^^- = oo, then Z^(p; z) is an entire function in both ar- 
guments z and p. 

(2) If limu^oo = b > 0, then I^(p; z) is analytic in the disc \z\ < ab and for any fixed 
p = 0, — 1, . . ., 2wi(/o; z), as a polynomial, is an entire function. 

For any p ^ 0, — 1, . . ., I^(p; z) admits an analytical continuation in the entire complex 
plane cut along the positive real axis given by 

1 f io ° ( z\ s 

(3.7) 2^(p;z) = — / a s (tp;a)T(s + p)T(-s) [ ds, | arg(-z)| < vr, 

2nrT{p) J_ ioo V a J 

where the contour is indented to ensure that all poles (resp. nonnegative poles) ofF(p+s) 
(resp. F(—s)) lie to the left (resp. right) of the intended imaginary axis and for any 
9te(s) > -1 

<p(a(k + s + l)) 



a s ((p;a) = Jj 



k=l 

with ip(s) = b — v(s) and v(s) = J °° e~ sr v{— oo, —r)dr is the Laplace transform of the 
tail of the Levy measure v. Moreover, for any p £ C, I^(p;z) admits, in the half-plane 
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VKt{z) < ff, the following power series representation 

(3-8) = (i - ^VtM-nW-^ (j^J • 

Finally, for any fixed «Re(z) < f , 2^(p;.z) is an entire function in the argument p. 

Proof. The first item and the first part of the second one follow from the discussion preceding 
the Proposition. Thus, let us assume that lim^oo^ = b > and p ^ 0, -1, . . . bmce 
ip S ip is well-defined and analytic in the positive right half-plane and vp(u) > for any 

u > 0. Our first aim is to extend the coefficients a n (ip; a) to a function of the complex variable. 
For the paths of the Levy process £ is of bounded variation, its Laplace exponent ip has the 
following form, see [31 Sec. VII. 3], 

ip(u) = u(b — v(u)). 

Thus, for any n > 0, we have 

a n (V;a) = TT7 — * a n (ip;a) 
T(n + l)a n 

with a n {tp;a)~ l = rifc=i c /'( a ^) an d ao(tp;a) = 1. It is plain that the mapping v is analytic in 
Fq = {s G C; 9\e(s) > 0} and v(u) is decreasing on M + with < u(0) < ab since ifi £ CK, + . 
Then, we may write 

a s (ip;a) = 1 s a s (f,a) 

1 (s + ljar 

1 ~ y(q(fc + S + l)) 

~~ r(s + l)a s A_l <p(afc) 

where the infinite product is easily seen to be absolutely convergent for any SHe(s) > by taking 
the logarithm and noting that \v(s)\ < v(9\e(s)). Moreover, a s (ip;a) satisfies the functional 
equation 

a 8+1 (tp;a) = —— — — -a s (tp;a) 
<p(a{s + 1)) 

which shows that a s (tp;a) is analytic in the half-plane F—\ = {s G C; 9te(s) > —1}. Con- 
sequently, a s ((p;a) is bounded on any closed subset of Then, we set G(s) = T(s + 
p)r(— s)a s (cp; a) and define 

where the integral is take in a clockwise direction round the contour consisting of a large 
semi-circle, of center the origin and radius R, lying to the right of the imaginary axis. This 
contour is intended to ensure that all poles (resp. nonnegative poles) of T(p + s) (resp. T(—s)) 
lie to the left (resp. right) of the intended imaginary axis. This contour is always possible since 
we have assumed that p / 0, — 1, . . . We can split 3£ R up into two integrals, '^% iR along the 
imaginary axis and, writing s = Re 1 ®, 

3 €r = — - / G(Re ie ) (--) Re i8 d6. 



2m J_ n / 2 v « 
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Recalling the following well-known asymptotic formulae, see e.g. \31[ Section 2.4], as |s| — > oo, 



T(-s) 



27re" 



Re ie R Re l9 +p- \ p i8(Re ie +p- §) 



2 e 



|f| <7T, 



e -7rR| sin 0| e Ke w ft-Re ie - 5 e id{-Re ie -i) 



< 7T 



and 



a 



~ | a 2| Rcose e _ - Rsin6 ' ar s(~ 2; ) 



we deduce that as |s| — ► oo 
(3.9) 



G(,)(-- V ^ 



[ sin0|(7r+arg(— z)) 

1 1 I g--«| sin0|(7r— arg(-z)) 



< 6» < vr/2 
-vr/2 < < 0, 



where a is a positive constant. On the one hand, along the path 2lj# we have = ±^ and thus 

as \z\ — > oo 



G(a) 



a 



~ a 



^(pJ-igif^Cp) 



e -i?(vr+arg(-2)) Q = ^ 

e -R(7r-arg(-z)) £ _ _ 7r /2. 



For the integral (|3.T[) to converge absolutely therefore requires that | arg(— z)\ < n. On the other 
hand, the asymptotic estimate (13.9H gives that as R — > oo 



- ► if \z\ < 1 and | arg(— z)\ < ir. 



Thus, as R — > oo, 

1 r io ° / z \ s 

55fL°w(-5)*- 

Finally, evaluating 3£ fl by the Cauchy integral theorem and letting i? — > oo , we get 

— / G(s) ( J ds = ——'^2a n (tp;a)T(p + n)z n , \z\ < 1 and | arg(-z)| < ir. 

■£nr J -ioo \ a/ l \p) 

Therefore, the integral (|3.7f) offers an analytical continuation of the mapping z 2(p; z) in 
the entire complex plane cut along the positive real axis. Moreover, we deduce from such an 
analytical extension that the power series (13. 6p has an unique singularity on the circle \z\ = ab 
located at the point z = ab > 0. Now, following a device developed for hyper geometric series, 
see e.g. [29], we introduce the function Ti defined for some a £ C by 



Hf ta (p;z) = (1- z) p; 



Note that 
(3.10) 



aabz 
z - 1 



Thus, denoting by (b n ) n >o the coefficients of the power series Ti.^ t a(p', z), we have &o = «o and by 
means of residues calculus, with £ a circle around of small radius and with positive orientation, 

li 



we have for n > 1 



b n = JL / Z^fU 

27TZ /* Z n+1 



oak \ <iz 



<- 1 )"i/,' 1 -^(^) 



Thus, one gets 



n n ~ ( actbz \ ^ , 



-n;aab) — —z . 



From Weierstrass's double series theorem, the above identity is true if jar] < 1 _A , . Moreover, 

the function on the left-hand side has a singularity at z = 1 and z = Thus, the series on 
the right-hand side is convergent if \z\ < 1 and \z(l — a)\ < 1. By choosing a = 1, we conclude 
by observing that the series on the right-hand side of (|3.10p is convergent for SHc(z) < □ 

Next, we introduce the following notation, for any z£C, 

O^z) = l^z) 

oo 

= ^(-ira n (V;a)z n , 

n=0 

and, similarly, for any p G C, 

0^{p;z) = l^(p;e m z), |arg(z)|<vr. 
Finally, in order to present our next result in a compact form, when HO holds, we write 

(3.11) j=\+® if * >0 ' 

I 8 otherwise, 

and, for any u > 0, 

\ipe\u) otherwise. 

We are now ready to state the main result of this part. 

Theorem 3.2. Assume that HO holds and recall that S(tx' a ) = P(xt~ & ) = Q X (T > t), x, t > 0. 
Then, there exists a constant C 7 > such that 

(3.13) 0^{ocy\x) ~ — — as x —> oo 
and 

(3.14) s(t) = c 1 r^o^(a r ,r 1 ) : t>Q. 

Moreover, the law of Tq under Qi is absolutely continuous with a density, denoted by s, given 
by 

s(t) = a 7 C 7 t- fi T- 1 0^(l + a7;t- 1 ), t > 0. 
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More generally, for any integer m, writing = J^-s, we have, for any t > 0, 

sM(t) = ( _ ir r(m + 1 + 07) c ^-x- mo ^ {m + 1 + ~ 7 . t -i y 

We postpone the proof of the Theorem until the Section 13,11 and offer instead a few conse- 
quences of this result. First, we deduce readily from the above Theorem the following asymptotic 
behaviors for large values. 

Corollary 3.3. With the notation used and introduced in Theorem \3. 21 we have 

(3.15) S(t) ~ C 7 t- 57 as t -> oo 
and, for any m = 0, 1 . . ., 

(3.16) S {m) (t) ~ (-l) m C £(^+l + ^Z) t -a7-l-m fls t^OO. 

T(ary) 

In the case lim^^oo = oo, by easily checking that the mapping cc h- > I^(p; x) is increasing 
on M + , for any ^He(p) > 0, one deduces, again from the above Theorem the following. 

Corollary 3.4. Let us assume that lim u ^oo = oo. With the notation used and introduced 
in Theorem \3.SX the entire function z \— > (dry, z) has no real zeros. 

Next, we observe that the function P £ B(A), hence P(xt~ a ), x,t > 0, is solution to the 
initial boundary value problem (12. 5ft . We do not aim to pursue in this direction, nevertheless 
let us point out that the uniqueness of the solution can be proved by invoking an argument 
involving the Martin boundary associated to the process X. Further details are provided in 
[39]. Next, let us recall that in [37], the expression of the Laplace transform of Tq, in the case 
E[£i] < 0, q = and 9 < a, is given, for any r, x > 0, as follows 

(3.17) E x [e~ rT °] = M^ e {rx a ) 
where 

N^ e (r)=l^{r)-C{e)r & % e (r) 
and the positive constant C(8) is characterized by 

(r) ~ C(0)r & % g (r) as r — oo. 

Thus, Theorem 13.21 offers a formula for the inversion of this Laplace transform. More specifically, 
one has 

/•oo 

M^ fi {rx a ) = x~ a / e- rt s(tx- a )dt. 



From the expression of s, it appears that a term- by-term inversion does not seem trivial. 

Moreover, as we have excluded the case when — £ is a subordinator, it is not difficult to 
verify that the Levy processes we consider in this Section satisfy the conditions listed below 
(|2.10p . Hence, from (|2.10p and (|3.15p . we get that the asymptotic behavior f|3. 13j) offers a 
characterization of the Kesten's constant. In what follows, we provide some representations of 
this normalizing constant in terms of the Laplace exponent. To this end, let us observe, in the 
case liniu^oo = oo, as < ip' 1 (0 + ) < oo, we have, for any u > 

ip^au) = bau+-{auf+ I (e aur - 1 - aur)e~< r v(dr) 
= {au) 2 (p 1 {au) 
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where b = b + 0-7 + f_ (e 7r — I{[t|<i>) fv{dr) and 

l /*oo r—r /*— 5 

</? 7 (cra) = — + - + / e- Qur / / e 7 V(cfo)dsdr. 
au 2 J J -00 J -co 



Thus, as above, one may define the function 



a s (ip 7 ;a) = j a s (y 7 ;a 

crl z (s + 1 



(- + 1) 

1 <fj(a(k + s + 1)) 

a 2 r 2 (s + l) AJL ^ 7 ( a fc) 

and observe the identity 

a s+ i((^ 7 ;a) = _ 1 a s (y 7 ;a) 
p 7 (a(a + l)) 

with ao((^ 7 ; a) = 1. Hence a s ((^ 7 ; a) is a meromorphic function in F_ 7 = {s G C; 9te(s) > — 7— 1} 
with simple poles at the points Sk = —k — 1 for k = 0, 1, . . . and Sk > —7 — 1. 

Proposition 3.5. J/lim^K, = 6 i/ien 

C 7 = a Q7 a_a 7 ((^ 7 ; a). 

Otherwise, we have, writing tp(au) = auip(au), 

y 7 (0+) i/a 7 = l 

C 7 = i a n ^ 7 (0+) (flLi "7 = n + 1, n = 1, 2 . . . 

> r( "_ a ~ 7 ) a-a 7 (^7;o0 otherwise. 

Finally, if q = and < # < a. Then, the constant Cq is such that 

1^ (r) ~ T(l - a9)C e r &e l^ e (r) as r -► 00. 

Proof. Let us start by pointing out that it is not difficult to check that we have, in all cases, 
C 7 > 0. Moreover, let us first assume that lim^^oo = b. From Proposition 13.14 we have 



1 f io ° ( z\ s 

0^(ar,z) = 2^(57) J a s(f-y;ct)T(s + aj)T(-s) ds, |arg(V)| 



< 7T. 



Hence, upon displacement of the path to the left in order to include the first pole of T(s + cry) 
we obtain, from Theorem 13.21 and a residue computation, that 



0^7 ("7! z ) = a Q7 a_a 7 (v? 7 ; a)z " 7 + o(z" 



■»7\ 



which gives the characterization of C 7 in this case. 

For the other case, one may follow a similar line of reasoning than in the proof of Proposition 
13. II to obtain, writing G(s) = "r^+'^ lXs + ory)r(— s), the following identity 

1 f ioc - ( z \ s 

^ ( « 7;z)= 2^^y_ 400 G(s) M ds 
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which is now valid in the sector | arg(z)| < tt/2. As above, after a displacement of the path to 
the left in order to include the first pole of T(s + a^) we obtain, from Theorem 13.21 and a residue 
computation, that 

([a-y] _ \ 

g »fjf^ {^y + a-^ow (±) ■ j + ^ 

where the sum is if [cry], the integer part of cry, is lower than 1. Since a s (ip^,a) has a simple 
pole at j = —1, . . . , —[cry], the terms in the sum vanishes. Hence, if is not an integer G(s) 
has a simple pole at — cry and the expression of C 7 follows readily in this case. If cry = n + 1, 
then G(s) has a double pole at — (n + 1) and using the recurrence relations of both the gamma 
function and a s ((^ 7 ;a), we deduce that 

Res s= _ (n+1) G(s) = s Jim ^ ((s + n + l) 2 G(s)) 




1 ± m _ I a " 2 II ^7( a ( s + k))^(a(s + n + l))a s+n+ i(^ 



fc=l / 
n 

= a" n - 2 r(n + 1)^(0+) 

fe=l 

and the result follows. The second part of the Proposition is proved as follows. Let us write 
F(r) = Ei [1 - e- rT °] . Then, from (IBTTj) . one deduces easily that 

F{r) ~ C(0)r 59 as r -» 0, 
which is equivalent, according to Bingham et al. [91 Corollary 8.1.7], to 

T(l — a0) 

which completes the proof. □ 

Note that the constant Cg also appears in the study of the rates of convergence of transient 
diffusions in spectrally negative Levy potentials. Indeed, Singh [421 Theorem 1] characterized 
the rates of convergence of these diffusions hitting times in terms of scaled one-sided stable 
distributions. It is easily verified that his scaling factors can be expressed in terms of Cg. Thus, 
Theorem 13.21 and Proposition 13.51 offer some explicit representations of these factors in terms of 
the Laplace exponent of the Levy potential. 

3.1. Proof of Theorem 13.21 According to Theorem 12. H our aim is to derive an expression for 
the function P(x) = 1 — Q x (Ho < oo), x > 0. We split the proof into two main stages. First, 
we show how to construct some time-space invariant functions for the semigroup of X when X 
does not reach a.s., i.e. X is associated via the Lamperti mapping with a ip £ £/C+. Then, we 
develop some devices for relating the previous result to the construction of the invariant function 
for the Ornstein-Uhlenbeck process defined in (|2. If) . To achieve this, we will consider separately 
the cases q > 0, and, q = and E[£i] < 0. 

We use the notation introduced in the Theorem 13.21 and take first X with underlying Levy 
exponent ipj. We denote its law (resp. its expectation operator) by (resp. iW). We 

recall that V7 6 C1C+ and hence the condition HO does not hold. Next, we simply write 
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Q(t) = (Qj ) t >o for the semigroup of X, i.e. for any bounded Borelian function g and t,x > 0, 
one has 

Q/ fl (x) = EW[ fl (X i )]. 

From [6], we have that Q 7 is a Feller semigroup. Then, we say, for any r£K, that a function / 
is r-invariant for Q" ( if 

e- rt Q]l{x) = I(x), x>0. 
We start with the following Lemma which is a slight generalization of [371 Theorem 1]. 

Lemma 3.6. For any x > 0, the mapping x t— > (x a ) zs an —1-eigenf unction for X, i.e., 
with the obvious notation, we have 

L^O lPj (x a ) + O^(x a ) = 0. 

Hence, for any r > 0, the mapping x i— ► 0^(rx a ) is —r-invariant for Q 1 . 

Proof. First, recalling that the mapping z \—* 0^ i {z) is an entire function, it is plain, from 
Proposition that the function 0^(x a ) £ D(L). Next, note that, for any f3 > 0, x 13 £ B(LT) 
and 

(3.18) LV 3 = x p - a il> y {P). 

Then, for any positive integer N, writing OS (x a ) = Yln=o(~ l)"' a n(V'7; a ) x ° Ln an d using suc- 
cessively the linearity of the operator L and (|3.18p . we get, for any x > 0, 



N 



V'Ol(x a ) = J> n (V 7 ;«)(-l) n L 



7 x ™ 



n=l 
2V-1 

= -Y^a n {^-a){-l) n x an . 

n=0 

The series being absolutely continuous, the right-hand side of the previous identity convergences 
as N — » oo. Hence, by dominated convergence, we get 

VO^{x a ) = -0^(x a ), x>0, 

which completes the proof of the first statement. Then, on the one hand, observe from the 
self-similarity property of X, that, for any r > 0, the mapping x t— > f{x) = 0^(rx a ) is an 
— r-eigenfunction. On the other hand, we have, for any r, t, x > 0, 

(3.19) j/ tQltf{x) = rert #/(*) + e rt j t QV(x)- 

Since Q 1 is a Feller semigroup, Qq/(x) = f{x) and 

j t Qlf{x) = QlVf(x). 

Thus, integrating f)3. 19|) yields 

e rt Qlf(x)-f(x) = r [ t e rs Q]f(x)ds+ f e rs Q^ f{x)ds. 
Jo Jo 

The proof is completed by recalling that L 7 /(x) = —rf(x). □ 
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Following a device developed by the author in [34J, we now show how to construct some 
specific time-space invariant functions for Cp in terms of the r-invariant function of X. We 
emphasize that this methodology may be applied to any self-similar Markov processes, that 
is also for the ones having two-sided jumps. Indeed, the proof of the following result, which 
generalizes a device suggested by Shepp |41j in the case of the Brownian motion, relies merely 
on the self-similarity property of X. Finally, we write, for any AgI, 



(3.20) 



+00 



if A > 0, 
otherwise. 



We are now ready to state the following result which is a generalization of [34, Theorem 1]. 

Lemma 3.7. Assume that, for any r < 0, the mapping x 1— > I{rx a ) is r-invariant for Q 7 and 
that I G L l (g(dy)). Then, writing \ = a ^ an d recalling that G K is a Gamma random variable 
with parameter k > 0, independent of Too, the mapping x \— > I R {x) defined by 

(3.21) I K (x) = r(«)E [I ( X G K x a )) , x > 0, 
satisfies the identity, for any < t < ( x , 

(3.22) (1 + X ty K Ql (d {1+xt) -aI K ) (x) = I K (x), x > 0, 
where d stands for the dilatation operator, i.e. d c f(x) = f(cx). 

Proof. Let us assume that / G L l (g{dy)). Then, for any x > and < t < £ A , we have, from 
the definition of I K , that 



(l + X t)- K Q/(d (1+X i)- a / K )(x) = (l + ^r K lEi 7) J°° l{xU^ + Xt)- 1 X?)e^y 



dy 



Performing the change of variable z = y(l + \t) 1 an d using successively Fubini theorem and 
the invariance property of I, we get 



[l+Xt)- K Ql (d( 1+xt )-a/«) (x) 



El 7 ) 



e- zxt I{ X zX?)e- z z K - l dz 



I { X zx a ) e~ z Z" 1 



dz 



Ux) 



which completes the proof. 



□ 



Next, we introduce the stopping times, and defined, for any A £ R and a > 0, by 

T a (A) = inf{s>0; X s = a{l + a\s)«} A ( x , 
ffW = inf{ S >0; e' Xs X Tx(s) = a} 
l)/a\. Writing (a)+ = max(a, 0), we have 



where t^(s) = (e aXs 
(3.23) 



and, in particular, for a = 0, since = Tq A ( x , we obtain 

e X«H^ (J (1 + aATo) + 
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In particular, we have H = Hq & ' with Hq defined in (|2.2|) . Next, for any a > 0, we set 

K +(a) = inf{K G M+; C^(k; a a ) = 0} 

K -(a) = mf{ K eR + ;l i ,(- K ;a a )=0}, 

with the usual convention that inf{0} = oo, for the smallest positive real zero of the function 
0^p(.;a a ) (resp. I^(—.;a a )). We are now ready to state the following. 

Corollary 3.8. Let < x < a and recall that \ = aX. 

(1) Then, for any A < and any p G C with 5He(p) < n^(\x\ a a), we have 

L J G^{p; \x\a a ) 

Consequently, for any real k such that k < (a), the mapping x i— > is 
positive on M + . 

(%) Similarly, for any A > and any p G C wif/i 9te(p) > {x a °)> we have 

z^(p; xx a ) 



(i + xtW)-" 



where we assume that Xa Q <C b if lim w _ 



i/> 7 (u) 



b. 



Proof. Let us first deal with the case A < and by the self-similarity property of X, we set, 
without loss of generality, \ = ~ 1- Moreover, since X is a Feller process on [0, oo), we can start 
by fixing x = and a > 0. Then, recalling that 0^ 7 (O,a a ) = 1, we observe that C?^, (/c;a a ) 
is positive for any < k < (a) reals. The existence of such an interval follows from the 
fact that the zeros of a non constant holomorphic function are isolated. Thus, by means of the 
Dynkin formula combined with Lemmae 13.61 and 13.71 we deduce, for any < k < k \ (a), that 



(3.24) 



E, 



(7) 

o 



1 



Next, we recall, from the identity (|3.23p . that 

e ^(J(l_ r (A))- K . 

bmce is a positive random variable, as a Laplace transform, the left-hand side on the 

identity (|3.24j) . is analytic in the half-plane {p G C; !SKe(p) < rej (a)} and positive on E + , 
see e.g. |44} Chap. II]. Then, let us assume that there exists a complex number p(a) in the 
strip < SRe(p(a)) < (a) such that G^{p{a);a a ) = 0. However, as the left-hand side of 
(|3.24p is analytical with respect to the argument k in this strip, we deduce, by the principle of 
analytical continuation, that this is not possible. Moreover, we get that (p; a a ) has no zeros 
on {p G C : 9te(p) < (a)} and is positive on {k G R : n < (a)}. Finally, let us consider 

a real number a± such that < a\ < a. Clearly, Q^-a.s. (1 — T^' )^ K < (1 — ?a for any 

< k < (cl) A (oi). Then, we deduce from (|3.24|) . for any < k < (a) A (ai), that 

1 1 



< 



< 



G^(K;a c 



IN 



Thus, it is not difficult to see that (g&i) > (a). Therefore, since (x) > «^ (a), for any 
< x < a, the strong Markov property and the absence of positive jumps of X complete the 
proof in the case A < 0. The case A > is proved by following a similar reasoning. We point out 



that in [34} Corollary 3.2], the restriction Xa a < b if lim u _ 



few 



b has been omitted. 



□ 



We are now able to proceed to the proof of Theorem 13.21 which we now split into two parts: 
the case when To occurs continuously, i.e. q = and E[£{\ < and the case when To occurs by 
a jump, i.e. q > 0. 

3.1.1. Continuous killing. Here, we assume that q = and E[£{\ < 0. Thus, there exist 9 > 
such that ip(9) = and 7 = 8. We recall that ipe(u) = ip(9 + u) with ip' e (0 + ) > and thus 
4)0 € CK+. 



Lemma 3.9. (1) Let A < 0. Then, K+(|x|"a) 
Kf(\x\ a a) and < x < a, we have 



(i + xtW)-i {t (, )<ToACA} 



(3.25) 



In particular, for any < x < a, we have 



< To A C ; 



(2) Similarly, for any A > and < x < a, 



(3.26) 



Tj A ) < To 



a9 > and for any k < 

x 9 Q^ g (K + a9;xx a ) 
a 9 O^ e {K + a0;xa a ) ' 

x 6 O i , e (a6;\x\x a ) 
a 9 0^ e {a9;\ X \a a )' 

x 9 l^ e {a9]xx a ) 
a 9 I^ g {a9- X a a ) 



b. 



where we assume again that \a a < b if lim u - 

Remark 3.10. From f)3.25j) and (|3.26p . we observe, for any A < and < x < a, the following 
identity 

Q x < To A C a1 



T^l a < T 



,(A) 



It would be interesting to prove such a formula directly from the definition of T„ 

Proof. Let us first consider A < and as above we fix x = -1 We start by using the fact that 
the function x 1— > x~ 9 is excessive for , see e.g. [40] . More precisely, one has, for any t > 
and for any F a J^-measurable and bounded random variable, 

Ef [F] 



X°F,t<T 



x > 0. 



Note that this relation also holds for any .Foo-stopping time. Moreover, proceeding as in the proof 
of Corollary 13.81 one gets that the Mellin transform of the positive random variable (1 — Ta^) + 
is well defined for any real k such that k < 0. Thus, since X has no positive jumps, one obtains 
by means of both Corollary 13.81 and the optional stopping theorem, for any k < 0, 



E, 



{Ti A) <T } 



£_ E w [(i-tW); 

x 9 0</>o( K + a8;x a 
a 9 0^(k + d8;a a 



-(n+aO) 
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We deduce that K^~(|x| a a) > and the proof of the item (1) is completed by letting k — > 0. The 
last statement follows again by adapting the above arguments to the case A > 0. □ 

We are now ready to complete the proof of Theorem 13.21 in the case 7 = 0. Let us fix A = —a 
and observe, from (|3.2U|) . that (~~ a = 1, then one gets that 



ri~ fi ) < T n A 1 



}x [r-a(H a ) < T_ 5 (# ) A 1] 
he [Ha < H ] 



since the mapping t *—> r^a{t) is increasing and t ~(1) = 00. Thus, as X has no positive jumps, 
one deduces that 



lim 



Q x [H = 00] 
P(x). 



Hence, as we have learnt from Corollary 13.81 and Lemma 13 . 9 1 that the mapping x 1— ► O^ e (a0;x a ) 
is positive on M + , there exists a constant Cg > such that 



Finally, recalling that lim a 



P(x) = 1, we obtain that 
p(x) = C e x e O^(a0;x a ). 



Hence, we deduce the expression of S by recalling that S(t) = P(t~ a ). Finally, the series 
O^ (a8; x a ) being absolutely continuous, the expression of the density s is obtained by differ- 
entiating terms by terms. Indeed, one has 



»(t) 



Cat 



-&6-1 



_. 00 



T(&9) 



l a n (V0(a# + n)Y(aO + n)i~ 



ra=0 



r(d# + i; 



The expression of the successive derivatives are obtained by an induction argument. 



3.1.2. Killed by a jump. Throughout this part, we assume that £ is a spectrally negative Levy 
process killed at some independent exponential time of parameter q > 0. Recall that, for any 
u > 0, tp(u) = ip(u) — q, <p is such that ip o <f>(u) = u and the mapping ip<p( q )(u) = ip(u + 4>(q)) € 

OC+ as V4( 9) (0+) = > 0. 
Lemma 3.11. Let A < 0. Then, ^ q) {\x\ & o) 



kZ (|y| a a) and < x < a, we have 



fii (|x| a a) — > and for any k < 

*V(<z) 



:i+xtW)^i {T (x )<To} 



a0(9) (K + d0( g );x« a )' 
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In particular, 



zf> < To A C ; 



Similarly, for any A > and < x < a, we have 



Q^(g^( g ); X x a ) 
j*(9) 0,/, (30(g) ; X a a ) 



where we assume that \a a < b if lim u _ i . 00 -^^-H- = 5. 

Remark 3.12. From the Lemma, we also observe, for any A < and < x < a, the following 
identity 



T a W < T A 



T->(W) / rp 



Proof. Let us observe from the Lamperti mapping (jl.ip that the semigroup, (Qt)t>o of X, is 
given, for a function / positive and measurable on R + , by 

Q t f(x) = 4«> [e-^/(X t )] , t > 0, x > 0, 

where stands for the expectation operator associated to the law of X with underlying 
Laplace exponent ip. This identity could also be derived from the Feynman-Kac formula. Thus, 
for any ^-"oo-stopping time T, one has 

E x [f(X T )]=E^ [e-* A Tf{X T )]. 

Moreover, as £ has independent increments, it is plain that the process (e~ qt+ ^ q ^ t )t>o is F^- 
martingale, where P^) stands for the law of the Levy process with Laplace exponent ip. By time 
change, one deduces that the process )t>o is a Qi-martingale. Thus, one can define 

a new probability measure, which we denote by <Q>(^w), as follows, for any t > and for any F 
a J-"t-measurable and bounded random variable, 



E(f(«)) [F] 



X, 



<t>(q)-qAt 



, x > 0. 



It is easily seen that the underlying Laplace exponent of X, under Q\fw) ^ is ipM q )- Hence, one 
gets by the absence of positive jumps for X and an application of the optional stopping theorem, 
that, for any < x < a and k < 0, 



E. r 



:i + xrW); K i {T (, )<To} 



E<«> fe- 9 V)(l + x TW);% rW 



{Tti A, <T }_ 
(A)y-(«+5^(g)) 



a/ 

xn*) c Vrt, ) (* + 5 ^(?);lx|3 a ) 



(D 



^ w (K + ^(9);lxl« a ) 



where the last line follows from Corollary 13.81 since ^^^(0 + ) > 0. The first statement is proved. 
The second one is obtained by means of the same devices. □ 

The proof of the Theorem is completed by following a line of reasoning similar to the previous 
case. 
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4. Some illustrative examples 



4.1. The Bessel processes. We consider £ to be a 2-scaled Brownian motion with drift 2b 6 R 
and killed at some independent exponential time of parameter q > 0, i.e. tp(u) = 2u 2 + 26n — q 
and 20(g) = ^/2q + b 2 — b. Note that 4><j>(q)( u ) = 2« 2 + (26 + 0(g))u. Its associated self-similar 

process X is well known to be a Bessel process of index b killed at a rate q X~ 2 ds. Moreover, 
we obtain, setting g = b + 20(g), 

" ' n=0 

= $(p, e +l;-x/2) 

where $ stands for the confluent hyper geometric function. We refer to Lebedev |25|, Section 9] 
for useful properties of this function. Next, using the following asymptotic 

$ (p, 8 + 1; -x) ~ r ( g + 1 ) x -p as x ^ 
r(e + 1 - p) 

we get that C^( g ) = ^i^v~pf) • Thus, we obtain, recalling that, for any g > 0, £ — 0(g) = 
6 + 0(g) >0, 

= *M n X)r(g + f ^- 1 * ( l + ^ b + 1; -Car 1 ) 

2^)r(0(g)) 7 V ; 

which is the expression [45 \ (5. a) p. 105]. Considering now the case q = and b < 0, we obtain 
readily that 6 = —b and 

a «(*) = ^^(i-M-fc;-^)- 1 ) 



2 6 



r(-6) 

dentity 

a Gamma random variable of parameter —b > 0. 



Hence, we deduce the well-known identity (To, Qi) = 2 g " wnere we recall that G_b stands for 



4.2. Law of the maximum of regular spectrally positive stable Levy processes. Here 
X is a regular spectrally negative a-stable, a 6 (1,2), Levy process killed upon entering into the 
negative half-line. Caballero and Chaumont [11 J characterized the characteristic triplet of the 
underlying Levy process. In [35), the author obtained the expression for its Laplace exponent. 
Instead of using this expression, we follow an alternative route. Indeed, in [32], the author 
computed the unique increasing invariant function, say P\, of the Ornstein-Uhlenbeck process 
defined by 

U t = e~ t X^ ea \t_iy a \, t - °> 
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for A > 0. The function P\, is given, with C a constant to be determined and recalling that 
a = 1/a, by 



(4.1) 



Px(x) 



n=0 



(an + a) 



Cx 



a-l 



2*1 



;i,i),(i,i-5) 

(a, a) 



a\x L 



x > 0, 



where 2*1 stands for the Wright hypergeometric function. We mention that P\(x) is obtained 
from the identity (j3.5|) with the so-called g-scale function of a spectrally negative a-stable pro- 
cess which was computed by Bertoin [1] in terms of the Mittag-Lefher function. More precisely, 
we have 



with 



P x (x) = Climr(/e)E 1^ (G^(aX) a x 



1^ (x) = J2 



n=0 



r(a(n + l))' 



From Remark 13.121 we have P(x) = Pa(e ian x). Note that P(0) = and using the large 
asymptotic of the function 2*1, detailed can be found in [36], we get, as x — > 00, 

-1 



2*1 



(1,1), (1,1 -a) 
(a, a) 



sin(d7r) 



7T 



X 



Hence, by setting C = sin ^ t7r ) ; we obtain the required condition lim a ._ >00 P(po) = 1 and 
P(x) 



7T 



(1,1), (1,1-5) 
(a, a) 



Next, from the identity (|2.7|) . we find that 



max Z s > x 

0<s<l 



P(x) 



where Z is a spectrally positive stable process of index q. Thus, by differentiating, one gets the 
following expression for the density 

sin(avr)^_ 2 / (1, 1), (1, 1 - a) 

(a, a-l) 



Pip) 



7T 



which is the expression found by Bernyk et al. [2J Theorem 1]. 

4.3. The self-similar saw-tooth processes. Finally, we consider the so-called saw-tooth pro- 
cess introduced and deeply studied by Carmona et al. |15| . It is a self-similar positive Markov 
process of index a = 1 with underlying Levy process the sum of a drift of parameter 6 = 1 and 
the negative of a compound Poisson process of parameter (3 > whose jumps are exponentially 
distributed with parameter 5 + (3 — 1 > 0, i.e. 

u + 6-1 



ip(u) 



u 
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u > 0. 



Moreover, in [15], the authors show that 

m = \ [q - (S - 1) + 0(g)) , q > 0, 

where 4>(q) = \/(q — (5 — l)) 2 + 4(<5 + (3 — l)q. Let us proceed with the case q = 0. Note, for 
1 — /? < 5 < 1, that 7 = 1-5 and 

u + 1 - <5 

Thus 

/ I ,s T(n + 1 + /?)r(2 - 5) 

anWl-6, 1) = r(1 + mn + 1)r(n + 2 _ s y «0-l, 

and for |z| < 1 

j " r(p)r(l+/3)^ 1 j T(n + 2-5)n! z 

= 2 F 1 (p,l + P,2-S;-z) 

where 2F\(a,b; x) stands for the hypergeometric function, see Lebedev |25^ Section 9] for a 
detailed account on this function. Next, recalling the identity 

P , ,, ax ^ r(2 - 8)T(n + 1 — 5 — 0) 
2 F l( -n, 1 + 0,6;!) = r (2 - J + n)r(l - J - fl ' 

we recover from (|3.8p the well-known identity 

z 



2 F 1 (p,l + [3,2-5;z) = (l-z)-P 2 F 1 (p,l-5-p,5; 



z - 1 



which provides an analytical continuation of the hypergeometric function into the half plane 
me(z) < |. Finally, using the following asymptotic 

«fl(M + ft»-*--)~ rp:?-^+flj x "' 



one obtains 



Moreover, after some easy computations, one gets for 7 = 4>(q),q > 0, 

7 / \ u + <j)(g) 

^ {u) = U u + [3 + 5 + m-l 

Thus, proceeding as above, we obtain 

0^ {q) (p;z) = 2 Fi(p,p + 5 + ij>(q),l + <Kqy,-z) 

and 
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